Surface-guided computing to analyze subcellular morphology and membrane-associated signals in 3D

Signal transduction and cell function are governed by the spatiotemporal organization of membrane-associated molecules. Despite significant advances in visualizing molecular distributions by 3D light microscopy, cell biologists still have limited quantitative understanding of the processes implicated in the regulation of molecular signals at the whole cell scale. In particular, complex and transient cell surface morphologies challenge the complete sampling of cell geometry, membrane-associated molecular concentration and activity and the computing of meaningful parameters such as the cofluctuation between morphology and signals. Here, we introduce u-Unwrap3D, a framework to remap arbitrarily complex 3D cell surfaces and membrane-associated signals into equivalent lower dimensional representations. The mappings are bidirectional, allowing the application of image processing operations in the data representation best suited for the task and to subsequently present the results in any of the other representations, including the original 3D cell surface. Leveraging this surface-guided computing paradigm, we track segmented surface motifs in 2D to quantify the recruitment of Septin polymers by blebbing events; we quantify actin enrichment in peripheral ruffles; and we measure the speed of ruffle movement along topographically complex cell surfaces. Thus, u-Unwrap3D provides access to spatiotemporal analyses of cell biological parameters on unconstrained 3D surface geometries and signals.

Cell surfaces are extracted from binary segmented image volumes and stored as a mesh, a data structure described by a list of the Cartesian 3D vertex coordinates on the surface and a second list specifying how the individual vertices are connected into triangles or faces. Tracking the correspondence between two 3D surface meshes is an active area of research in computer graphics [25][26][27][28] and none of the methods have been adopted to cell imaging. A particular technical challenge that arises when adapting techniques from computer graphics with applications to cell biology is the non-convexity, irregularity and high curvature of surface protrusions on most cell shapes. Very few methods have been proposed to accurately follow such geometries over time and have largely been demonstrated on well-defined shapes such as human pose 29 or hands 30,31 . Generally, these methods track by matching meshes from consecutive timepoints. To match meshes, methods attempt to assign a unique signature per vertex or face to establish a matching between vertices and faces by minimizing a loss metric 32,33 . However, this approach is inherently sensitive to mesh quality, uniqueness of the signature, optimizer convergence and is difficult to generalize when tracking surfaces over many timepoints. Crucially, meshes segmented from two different timepoints have different numbers of vertices and faces and the lack of the exact same surface features poses ambiguity in matching. Alternatively, individual 3D image volumes may be first registered spatiotemporally before mesh extraction, as is done in neuroscience [34][35][36][37] . This approach is robust to mesh quality and sampling errors, however deformations must be small between timepoints. For example, long or thin cell surface structures such as lamellipodia and filopodia suffer voxel undersampling, limiting the registration to cell surfaces with largely globular features such as blebs 38 . This problem could be remedied by mapping the 3D surface to the unit sphere [39][40][41] . Indeed, in macroscopic imaging application, this procedure has enabled registration of complex geometrical features such as brain folds directly on the 3D sphere or in derivative 2D unwrapped images 27,[42][43][44] . Unfortunately, this strategy requires closed surface topologies with no holes (genus-0 surfaces), which is not generally guaranteed in live-cell microscopy. Alternatively, one can selectively segment surface motifs and track these in 3D whilst mapping surface-proximal molecular signal intensity 38,[45][46][47] . Like the 3D surface tracking, this weighted principal component analysis of surface variables, such as the local curvature (Extended Fig. 1f,g, Methods). If the input Cartesian 3D surface mesh ( , , ) is genus-0, the generation of a reference surface, ref ( , , ) may be skipped and 2D equiareal surface unwrapping realised directly (Fig. 1b, Suppl. Video 2). We note that an input genus-X Cartesian 3D surface ( , , ) is also directly unwrapped into 2D through steps 1-4 of u-Unwrap3D, but not in an equiareal manner (Extended Fig. 1h). In Step 5, the first part remaps the Cartesian 3D volume ( , , ) and associated signals ( ( , , )) into a topographic volume ( , , ) coordinate system that is normal to the reference surface, ref ( , , ). The second part establishes a bijective mapping of ( , , ) to ( , , ) coordinates, ref ( , ) , i.e. the ( , ) parameterized reference surface of ref ( , , ), is propagated in Cartesian 3D space in the surface normal direction at equidistant steps of voxels along the steepest gradient of the signed distance function, ∇Φ( , , ) , for a total of steps. Interpolation of the respective Cartesian volumetric signal intensities, ( ( , , ) at the ( , , ) coordinates indexed by ( , , ) generates the topographic 3D volume equivalents, ( ( , , )). Finally, in Step 6 the topographic 3D surface representation, ( , , ) of the input surface ( , , ) is obtained by surface meshing the topographic binary volume segmentation.
In summary, u-Unwrap3D generates bijective mappings of a given genus-X surface between 5 equivalent surface representations; Cartesian 3D, ( , , ), genus-0 reference 3D, ref ( , , ), the unit 3D sphere, 2 ( , , ) , topographic 3D, ( , , ) , and the 2D plane image, ( , ) (Fig. 2a), while simultaneously transforming Cartesian 3D to topographic 3D volumes (Fig. 2b). This was made possible by two crucial choices; the use of cMCF and voxelization to construct a genus-0 reference surface, ref ( , , ) to realise spherical parameterization (Step 1) and the implementation of an efficient numerical scheme to relax area distortion on the 3D sphere (Step 3). The former allows us to construct ref ( , , ) as a proxy of the genus-X ( , , ) surface mesh and to unwrap this 3D surface into one 2D ( , ) image, instead of requiring multiple 2D ( , ) images, which simplifies downstream analysis [56][57][58][59] . The latter ensures that the unwrapped 2D ( , ) image captures the salient surface features of ref ( , , ), and by extension the genus-X ( , , ) surface. Importantly, the bijectivity of the mappings guarantees that for any point on any of the surface or volume representations matching points exist on any of the other surfaces or volumes. Moreover, the bijectivity guarantees preservation of the point topology, i.e. a series of points ordered in clockwise fashion on one surface representation maps to a series of points ordered in the same way on any of the other surface representations and preserves the local neighbourhood relationships. As a result, we can apply mathematical operations defined in any one of the representations and map the results to any other. u-Unwrap3D thus supports the optimal spatiotemporal analysis of unconstrained surface geometries and associated signals.

Validation of u-Unwrap3D on diverse surface motifs
We validated the generality and performance of u-Unwrap3D by application to 66 single cell images acquired by high-resolution light sheet imaging 45,60 . The dataset span morphologically diverse cells with blebs, lamellipodia and filopodia. The cell surfaces were meshed with marching cubes and segmented within the u-Shape3D software 45 . Small errors in the initial segmentation and meshing process cause high-order genus surfaces with topological holes and handles, which cannot be unwrapped directly (Extended Fig. 2a). Holes can also generate non-watertight surface meshessurfaces that are not closed and have no clearly defined inside volume 48,61 possessing potentially complex internal volumetric structures that violate the assumptions of standard 3D mesh processing algorithms.
We first tested the number of input cell surfaces for which u-Unwrap3D could successfully run all steps 1-6 and compute all 5 of the representations as a measure of generality and robustness. Notably only 6/66 (11%) input cell surfaces were genus-0 and only 36/66 (55%) were watertight (Extended Fig 2b). In 63/66 cases, (>95%) we successfully ran all steps and obtained all representations (Extended Fig 2b). The three failures occurred in scenarios, in which the holes and handles remaining after the application of cMCF were still too large for the volume dilation to generate a genus-0 reference surface after remeshing (Fig. 1b, Step 1) (Extended Fig. 2c). In all successful cases, cMCF and binary voxelization under volume dilation generated genus-0 reference surfaces within a median of 10 iterations (Extended Fig 2c, c.f. lamellipodia). Fig. 2c shows extracted representations for challenging examples with blebs, lamellipodia and filopodia (Suppl. . We next tested the robustness and performance of the ref ( , , ) spherical parameterizations, (Fig. 1b, . Extended Fig. 3a confirms that the quasi-conformal spherical parameterization (Step 2) minimizes the conformal error to the ideal value of 1, with the largest error in cells with filopodia (1.016±0.013). We also verified the need to relax local area distortion. Whilst quasi-conformal spherical parameterization Ω 2 ( , , ) is equiareal for blebs (0.978 ± 0.037), the median area distortion showed that the surface fraction of lamellipodia was down to 0.432 and in filopodia to just 0.140 with respect to their original area fraction on the reference 3D, ref ( , , ) surface, let alone ( , , ). Our scheme for area distortion relaxation (Methods) produces a quasi-equiareal spherical parameterization (Step 3) in blebs (0.985±0.012), and successfully achieves the ideal value of 1 in lamellipodia (1.000±0.000) and filopodia(1.000±0.001), within a maximum median of 23 iterations for lamellipodia (Extended Fig.3b, Table i). We further tested the ability of our relaxation scheme to balance the trade-off between the two extremes of conformal to equiareal spherical parameterizations using different stopping criteria (Extended Fig.3b, Table ii-iv). The initial parameterization without any area-distortion relaxation (iteration 0) is by design conformal but also found to satisfy the most isometric parameterization (MIP) 62 . Running for Ω < a maximum of 50 iterations yields an equiareal parameterization for all motifs. At ≈ 1 2 Ω iterations the relaxed mesh jointly minimizes the summation ( + ln , Methods) of conformal ( ) and area distortion ( ) errors. At ≲ Ω iterations the relaxed mesh is the area-preserving MIP 63 . As expected, this latter parameterization does not fully minimize area distortion (blebs (0.997±0.011), lamellipodia (0.979±0.005) and filopodia (0.980±0.028)) but exhibits slightly lower conformal errors and consequently higher quality faces than a pure equiareal mapping.
Lastly, we tested how accurately topo ( , , ), which defines the topographic 3D mesh, ( , , ) (step 5) mapped back into Cartesian coordinates reconstructs the input surface, ( , , ). For all cells, ( , , ) was computed with a ( , ) image grid size of 1024x512 pixels. The aspect ratio, 2 x ( = 512) was chosen to preserve the ratio between the equatorial circumference and the length of the arc between north and south poles of a sphere. Compared to the input surface ( , , ), ( , , ) is lower genus and provides higher face quality (Extended Fig. 3c). We assessed the discrepancy between topo ( , , ) and ( , , ) using 4 metrics; Chamfer distance (CD), sliced Wasserstein distance ( 1 ) 64 , and differences in total surface area (Δ ) and volume (Δ ) (Extended Fig. 3d, Methods). Considering inevitable rasterization errors when mapping the floating-point precision 3D sphere Ω 2 ( , , ) to ( , ) defined on an integer ( , ) image grid, we measured low vertex position errors according to CD and 1 . Cells with lamellipodia had the lowest error (median CD=1.77 voxel, 1 =0.93 voxel) and cells with blebs were slightly worse (median CD=2.79 voxel, 1 =4.28 voxel), likely due to their intrinsically small height (small topographic ). As one would expect, cells with long, thin filopodia displayed the largest discrepancies (median CD=10.33 voxel, 1 =18.45 voxel). Correspondingly we measured a small Δ (+1.2%) and Δ (+7.9%) for cells with lamellipodia. Δ was larger for cells with blebs (-11.6%) and measured to be too large for filopodia (-55.3%) when compared to Δ differences measured after making ( , , ) watertight (+4.2% blebs, +3.4% filopodia). Visualization of exemplar cells show good geometric correspondence between topo ( , , ) and ( , , ) (Extended Fig. 3e). Salient surface features were largely captured, albeit smoothened and blurred in topo ( , , ) when local surface regions were underrepresented due to being distant relative to ref ( , , ) (Extended Fig. 3e, black triangles, 1 st row blebs and 4 th row lamellipodia). Most of the primary morphological features, namely the length and thickness of long, thin filopodia (Extended Fig. 3e, green triangles), except those located both densely together and distant relative to ref ( , , ) (Extended Fig. 3e, red triangles), were captured. u-Unwrap3D was able to capture both the complex lamellipodia folds and curved cell bodies to high accuracy (Extended Fig. 3e, row 2,4). Closer inspection of topo ( , , ) and ( , , ) in these cells traced a large Δ to meshing errors in the input surface ( , , ), which caused internal volumetric structures to be merged into the cell surface representation, and overestimation of total surface area. These errors affect the CD and 1 to lesser extent.
In summary, our results demonstrate that u-Unwrap3D is robust and applicable to process unconstrained geometries. For maximum resolution of high curvature surface features, a genus-0 reference surface ref ( , , ) proximal to the input surface ( , , ) is recommended with a large ( , ) grid size and small step sizes when propagating ref ( , ). However, these choices depend on the quality of the input surface mesh ( , , ), which depends on the robustness of cell segmentation in the face of noisy image raw data.

u-Unwrap3D enables unsupervised instance segmentation of subcellular surface motifs
The unbiased identification and segmentation of individual protrusive features in unconstrained 3D surface geometries is nontrivial. Cellular protrusions present complex morphological characteristics that are difficult to define descriptively. Even well-known morphological motifs exhibit significant heterogeneity and ambiguity. Not all blebs are spherical, lamellipodia are often plate-like with high curvature ridges but otherwise have no readily-defined shape prior, and filopodia, though long and thin, can sprout haphazardly from elevated 'stumps' or even off of each other. In areas of dynamic and dense protrusions, where does one protrusion start and another end? Consequently, most existing approaches focus on particular surface features of interest such as 'ridge' networks that can be segmented by designed imaging filters or through trained semantic segmentation, with morphological processing and parameter tuning 46,47,65,66 . With u-Shape3D we introduced a multi-class morphological motif detection by partitioning the 3D surface into convex patches and applying support vector machines trained with expert annotation to classify the patches into pre-specified motif types 45 . However, this approach cannot detect and segment all protrusions generally, only the limited motifs for which the supervised classifier has been trained on. Lastly, even after obtaining the segmentations, how do we systematically measure salient protrusion properties? For example, with respect to what reference surface should protrusion height be measured? Where is the protrusion width to be measured? How is the internal volume of a protrusion determined?
These segmentation and characterization problems can be significantly better defined in the topographic 3D surface representation ( , , ), which captures in one field-of-view all surface features protruding normally to ref ( , , ). As preserves the total Cartesian 3D curvilinear distance from ( , , ) to ref ( , , ) along the gradient of steepest descent we can formally define a 'protrusive' feature as having a -coordinate greater than that of a reference topographic surface, ref ( ref , , ), and measure the protrusion height as the difference, ℎ = − ref . For example, protrusive features could be specified as those having ℎ > ℎ ̅ , the mean height of ( , , ) relative to the planar topographic 3D cMCF surface  (Fig. 3a,b). Whereas this problem is difficult to frame in Cartesian 3D, in the topographic space the interpolation can be solved naturally by using asymmetric least squares (ALS) optimization with a Whittaker smoother 67,68 , where the asymmetric weights allow us to account for the heterogeneous protrusion height; and the desired level of surface smoothness can be incorporated as a regularization term (Extended Fig. 4b, Methods). To use ALS, we create a ( , ) -parameterized approximation of ( , , ) with ≃ ( , ) using straightforward image processing procedures (Extended Fig. 4b, Methods). By exploiting these properties of the topography space, we developed a general approach to segment any protrusion motif including blebs, lamellipodia and filopodia (Extended Fig. 4c,d) with minimal heuristic parameters to tune. Importantly, we did not need to design specialized image filters 65,66 , compute and cluster feature descriptors 46,[69][70][71] , or require data training [45][46][47] .
We demonstrate the segmentation of individual protrusion instances, capturing motifs identified by uShape3D, but without the need for training annotations. By construction, in topographic 3D ( , , ) space all surface protrusions are oriented upwards with increasing . Moreover, the tops of protrusions are individually separated as local regions of high topographic mean curvature, which we identify by thresholding and applying connected component labelling in the topographic volume, ( , , ) (Fig. 3c). Mapping the segmented regions back onto the topographic mesh ( , , ), we diffuse these initial 'seed' labels across the surface using a combined geodesic distance and dihedral angle affinity matrix to naturally segment the 'stem' of the individual protrusions (Methods). The dihedral angle measures the discontinuity in local mean curvature. It incorporates the prior intuition that the boundaries of a label should expand faster on local surfaces of homogeneous curvature such as that on a 'hill', compared to another label experiencing large curvature differences in its local surface region, such as in a valley between multiple 'hills'. The combined affinity matrix thus introduces a morphology-aware competition between segmentation labels and provides a biophysical rationale for defining which surface patches belong to individual 'seed' protrusions. Furthermore, the dihedral angle is large between a protrusion and the main cortical cell surface and thus serves as a soft stopping criteria for diffusion (Extended Fig. 4d) in addition to applying the binary protrusion segmentation from above. Lastly, we filter out protrusions that are too small and close any small holes using the Cartesian 3D surface area. The final segmentation result qualitatively and quantitatively agrees with that obtained by supervised u-Shape3D for lamellipodia ( Fig. 3d) but yields more contiguous labels and is less prone to oversegmentation (Fig. 3e). Importantly, this segmentation strategy is applicable, even when not all protrusions are equiareally represented in ( , , ), as shown by the segmentation of the majority of blebs and filopodia in exemplar cells (Extended Fig. 4d).
Individual segmented protrusions are genus-0 open-surface 3D submeshes that can be directly mapped to the 2D plane (Fig. 3f). This allows us to further refine the segmentation, for example, by detecting and splitting under-segmented blebs by a gradient watershed algorithm (Methods). Thanks to bijectivity, the refined segmentation labels can be transferred back onto the original surface mesh (Fig. 3f,g, c.f. before and after refine, black triangles). Both the before (adjusted normalized mutual information, NMI=0.57) and the after refinement (adjusted NMI=0.54) segmentations agree with u-Shape3D. However, u-Unwrap3D segmented blebs are more complete, with more blebs of larger surface area (150 blebs in total). In contrast, u-Shape3D over-segments small blebs (742 blebs) that were found to originate from erroneous meshing of internal structures in ( , , ) (see Extended Fig.3e). This example illustrates the potential pitfalls of identifying motifs from local surface patches only with potentially imperfect input 3D meshes. In contrast, with u-Unwrap3D any surfaces internal to the cell volume are readily removed when mapped into topography as these surfaces have -coordinates less than the reference surface, ref ( ref = smooth ( , ), , ).
Finally, the representation ( , , ) enables partitioning of the input cell volume into the sum of a reference volume representing the underlying cortical cell body and the unique volume occupied by individual protrusions (Fig. 3h). We do so by ( , )-parameterizing the reference cortical surface submesh, ref ( , , ) after removing all individual protrusion submeshes as a grayscale image such that the pixel intensity value at ( , ) is the respective -coordinate, and inpainting the missing -coordinates at ( , ) coordinates corresponding to the subtracted protrusions to generate a full reference binary volume (Extended Fig. 4e, Methods). For the protrusions, we first devised a marker-controlled lateral watershed depth propagation to diffuse the surface-based protrusion segmentation labels uniquely throughout the full topographic ( , , ) space slice-by-slice, top-to-bottom (Extended Fig. 4f). Then, individual protrusion volumes were generated by masking the propagated label volume with the reference binary volume (Extended Fig. 4g). We compared our topography-guided decomposition strategy to a fully Cartesian 3D mesh processing approach whereby individual protrusion submeshes were first closed by constructing a surface patch that minimized the local bending (or harmonic) energy 72 (Methods). The closed reference volume was then generated using all such patches to impute the residual holes in the raw reference surface. Fig. 3i panels i-ii show that the computed surface area (slope=0.90) and volume (slope=0.95) of individual protrusions are similar for both methods. However, the 3D mesh processing protrusions consistently under-measure larger protrusions. Crucially, the imputed reference volume appears artefactual. Where protrusions were located, the surface is overly smooth, and even involuted. These regions contrast starkly with non-imputed surface areas between protrusions, creating artificial 'peaks' and 'ridges' of high mean curvature (black arrows). In comparison, the topographyguided reference volume is mechanically more plausible.
In summary, the ability to map freely between topographic and Cartesian 3D surfaces and their respective volumetric representations enabled us to design simple and generalizable methods to detect and segment in an unsupervised fashion individual morphological motifs from unconstrained surface geometries. To bijectively map the topographically segmented surface protrusion labels onto the 2D plane, which is an optimal representation for tracking the segmented motifs, we developed a topographic cMCF for u-Unwrap3D (Fig. 3j, Suppl. Video 6, Methods).

u-Unwrap3D enables tracking of 3D subcellular surface motifs and molecular activity in 2D
A central goal of live cell imaging in 3D is to visualize the spatiotemporal relations between molecular activities and cell behaviors, including morphodynamic outputs. Progress has been made on software developments that allow unbiased and statistically meaningful analysis of cell morphology and molecular distributions 45,47,60,73,74 . However, to remain algorithmically and computationally tractable these pipelines have been restricted to quasi-static representations of dynamic processes. Surface-guided computing with u-Unwrap3D allows us now to remedy this limitation. Dynamic behaviours on complex 3D cell shapes, including their morphological and molecular signal activity changes, can be mapped to 2D representations where powerful analytical pipelines exist for spatiotemporally consistent tracking. Results can be statistically evaluated and, if of interest, be visualized in 3D by leveraging the bijective properties of u-Unwrap3D mappings. Fig. 4 and 5 demonstrate this capacity based on two examples.

Individual blebs dynamically recruit Septins to local surface regions during retraction
We first analysed potential relations between dynamic surface blebbing and the recruitment of Septins. Blebs are globular membrane protrusions of 1-2 μm diameter that are thought to extend in areas of localized membrane detachment from the actin cortex 75,76 . Intracellular pressure expands the budding blebs outward, followed by rapid assembly of a contractile actomyosin network that yields retraction. Cycles of protrusion and retraction have been described to last a few tens of seconds. Associated with the cycles are molecular activities both driving and responding to the morphological dynamics. One such process is the assembly of Septin protein polymers at sites of negative curvature (from a cell-external perspective) emerging at the bleb necks. Our previous work 38 has shown that disrupting the bleb cycle diminishes Septin assembly at the cell surface. Here we now exploit the ability of u-Unwrap3D to track individual bleb cycles and quantify Septin accumulation by remapping surface morphology and a fluorescent marker of Septins to an appropriate 2D representation. We acquired 3D volumes of SEPT6-GFP-expressing MV3 melanoma cells every 1.2s for 200 timepoints. As the cortical cell body exhibits little temporal variation and blebs protrude normally to the surface, the temporal mean cell surface, ̅ ( , , ) is a good proxy of the cell cortex. We apply u-Unwrap3D to ̅ ( , , ) to create a common static ( , , ) coordinate space for computing topographic 3D ( , , ) and 2D planar ( , ) representations for each timepoint. ̅ ( , , ) was computed by meshing the mean binary volume across individual binary voxelizations of the cell surface over all 200 timepoints (Methods). Note the construction of a common topographic space from a single mesh for a timelapse is computationally efficient but applicable only if cell shape changes lie within the Cartesian subvolume mapped by the topographic space. For large shape changes u-Unwrap3D should be applied to individual timepoints and spatiotemporal registration used to align surfaces to a common reference using ref ( , , ) or ( , ) representations. The segmentation tools discussed above were used to detect all bleb instances from ( , , ) and mapped to ( , , ) and ( , ) for each timepoint (Fig. 4a, Suppl. Video 7). Similarly, the computed mean curvature and the normalized Septin intensity surface signals ( ( , , )) from Cartesian 3D were mapped to topographic 3D ( ( , , )) and into 2D, ( ( , )), to enable simple bleb tracking and timeseries analysis.
To track blebs in 2D we computed the bounding box of individual bleb instances in every timepoint after appropriate image padding to account for spherical periodicity (Fig. 4b left). In ( , ) bleb dynamics can readily be followed by an established 2D multi-object bounding box tracker 77 with mean curvature optical flow-guided bipartite matching ( Fig. 4b middle, Methods). The bijective mapping allows us to map the resulting trajectories from ( , ) coordinates (individually colored) to ( , , ) coordinates to generate bleb tracks in Cartesian 3D ( Fig. 4b right, Suppl. Video 7). Due to the fast temporal acquisition rate, only a northern portion of the cell can be maintained in-focus. Again, u-Unwrap3D's bijectivity between 3D and 2D representation enabled us to easily map a manually annotated out-of-focus subvolume onto ( , , ) and into ( , ) to retain for analysis only the bleb tracks that remain within the in-focus surface regions (Fig. 4c). The distribution of individual bleb track lifetimes showed a peak at 14s and a long tail up to 240s, suggestive of a mixture of short-and long-lived blebs (Fig. 4d). The mean bleb lifetime of 27s corresponded well with a 30s periodicity given by the first peak of the temporal autocorrelation of mean curvature ( ( , , , )) in Cartesian 3D (Fig. 4e). This validates at the population level the accuracy of single bleb tracking after projection and segmentation of mean surface curvature in 2D.
The temporal autocorrelation curves of mean curvature and Septin intensity showed a high level of similarity, suggesting co-fluctuation of the two surface signals. Indeed, we had previously shown that surface regions of high Septin intensity with negative surface curvature for at least 30s display a correlation between negative curvature value and Septin intensity 38 . Whilst the majority of Septin pulses endured only one cycle of bleb formation and retraction, de novo formation of stable Septin structures appeared to be driven by several Septin pulses occurring in short succession. We thus hypothesized that these were formed by iterative blebdriven curvature generation events resulting in local levels of Septin oligomers surpassing a threshold necessary for inter-oligomer polymerization and enabling stabilization through formation of higher-order structures 38 . The ability to spatiotemporally track individual blebs enabled us now to quantitatively test this model. For the duration of each tracked bleb, we sampled within the 2D bounding box distortion-corrected timeseries of bleb surface area, on-/off-bleb surface mean curvature and Septin intensity (Methods). We used the 2D ( , ) bounding box to define a bleb's spatial area-of-influence and its surface area. The Cartesian 3D ( , , ) box area was taken as the bleb's 3D surface area in each tracked frame. Within the 2D bleb bounding box, 'on-bleb' is the largest spatially contiguous region of high mean curvature. The remainder area is 'off-bleb'. A single curvature threshold was computed by 3-class Otsu thresholding over all ( ( , , )) to define the regions of low/high mean curvature in ( , , ). Using the extracted timeseries we reconstructed the temporal profile of bleb area, on-/off-mean curvature and Septin intensity of a single mean bleb event in a window of 35s centered on the timepoint of maximum bleb area averaged over 545 single bleb events from 480 bleb tracks. We sample ±17.5s before and after the timepoint of maximum bleb area to exceed the 30s periodicity inferred from temporal autocorrelation by 5s and capture the full dynamics. We then broke the mean timeseries into four distinct temporal phases of bleb-mediated Septin recruitment, each ≈5s long (Fig.4f labels A-D). In phase A, the bleb begins to expand, increasing its surface area, accompanied by a sharp increase in on-bleb and a decrease in Septin intensity as the plasma membrane detaches from the actin cortex. The expansion also reduces off-bleb , causing a decrease in Septin intensity off-bleb, presumably due to disrupting Septin structures. In phase B, the bleb reaches maximum size and then begins to retract with decreasing surface area. Interestingly, unlike mean curvature, which begins to decrease before the maximum bleb area, the change in area is symmetrical, occurring ±2.5s relative to the time of maximum bleb area. Coincident with the bleb increasing to a maximum area, off-bleb decreases to a minimum and Septin intensity both on/off bleb stabilizes at a minimum. As the bleb retracts and the actin cytoskeleton reassembles, off-bleb increases and Septin intensity begins to increase both on/off bleb. In phase C (+2.5s to +8.0s after peak bleb area), the bleb area continues to decrease but at a slower rate. Unexpectedly, the off-bleb negative curvature plateaus at a value lower than its starting value (before phase A) instead of continuing to increase. Concurrently, Septin intensity undergoes the greatest rate of increase on both onbleb and off-bleb surfaces such that at the end of phase C, the Septin intensity is at the levels before phase A (c.f. -15s to -10s). In phase D (+8.0s to +14.0s), bleb area and on/off-bleb recover to pre-expansion levels while Septin intensity continues to increase before plateauing on both on/off-bleb surfaces. Beyond phase D, the next bleb cycle begins, with similar temporal characteristics to phase A. Altogether these results indicate that blebs generate optimal curvature dynamics during retraction to recruit Septin polymers to the surface regions around blebs. Moreover, the data support our model of Septins being recruited to negative curvature patches in a cyclic fashion, where each bleb formation-retraction drives the local accumulation of a few more oligomers until a threshold concentration is reached to trigger inter-oligomer polymerization of a stable Septin assembly. Notably, none of these observations could have been made without u-Unwrap3D.

Ruffles are driven by locally enriched filamentous actin and migrate actively on the cell surface
We also applied the capacity of u-Unwrap3D to examine putative relations between the dynamics of membrane ruffles, filamentous actin and surface actin retrograde flow. Membrane ruffles are thin, rapidlymoving, actin-rich protrusions 78,79 . They are thought to play a role in cell migration, and it has been proposed that ruffles arise as a consequence of inefficient adhesion in cellular lamellipodia 80 . Yet, it is unclear whether they have a specific function in migration or elsewhere, and the precise molecular and mechanical mechanisms of ruffle formation remain poorly understood. Although membrane ruffles have been the showcase for many of the recent advances in volumetric light sheet microscopy, there are no tools for quantitative assessment of membrane ruffling 78,81 . The thin, lamellar appearance make ruffles extremely challenging to segment in 3D live-cell images 46 . Moreover, ruffles are transient and they exhibit significant heterogeneity in nature and distribution across the cell surface. These characteristics make ruffles difficult if not impossible to directly track in Cartesian 3D.
There are many unanswered questions about cell membrane ruffles. Here, we focus on defining a membrane ruffle by a set of objective criteria and then on measuring their speed. We acquired 3D volumes of SU.86.86 pancreatic ductal adenocarcinoma cells plated on fibronectin-coated cover glasses and co-expressing Tractin-mEmerald and myristoylated CyOFP1. Tractin is a marker for actin 82 , while myristoylated CyOFP1 served as a diffuse cell membrane marker. We acquired images every 2.27s for 30 timepoints, and used u-Unwrap3D to map the Cartesian 3D ( , , ) surface of every timepoint into topographic 3D, ( , , ) and into 2D, ( , ) for tracking (Fig. 5a, Suppl. Video 8). The cortical shape change was relatively small but ruffles travel on the surface. Like for bleb tracking we constructed a common static topographic coordinate space ( , , ) for all time points. Here, we used the reference surface, ref ( , , , = 0) for this purpose (Extended Fig.5a). We then sought to track the spatial location of ruffles as individual 'ridge' objects. To do so, we projected ( , , ) onto the 2D plane, i.e. the ( , ) coordinate was the same as setting = 0, ( = 0, , ) instead of applying topographic cMCF. We computed the local enrichment of the actin signal as the ratio of Tractin-mEmerald to CyOFP1 intensity. Kymograph visualization of ( ( , , )) for mean curvature, and Tractin-mEmerald/CyOFP1 (TC) showed strong co-fluctuation of curvature and local actin intensity. It also highlighted the transient, 'ripple-like' nature of ruffles and their merging and dissipation with a delay between successive ruffles of ≈20s (Fig. 5b). To avoid complex image processing operations such as merging and splitting we applied optical flow-based region-of-interest (ROI) tracking 83,84 to the TC intensity to track simultaneously the protrusive membrane ruffling and retrograde surface actin flow (Fig. 5c, left). The ( , ) image grid size was 1025 x 512 pixels (not 1024 x 512 as ruffles required active contour cMCF to be used to generate ( , , ), see Methods) and the average dimensions of the tracked ( , ) ROI size was 23 x 23 pixels (approximately 0.06 x 0.11 μm pixel size in physical space). The resulting 2D ROI tracks exhibit unidirectional motion towards the cell center as shown by coloring directionality and remapping of ( , ) ROI tracks to polar ( , ) (Fig. 5c, middle, Extended Fig. 5b) and Cartesian 3D ( , , ) coordinates (Fig. 5c, right) (Suppl. Video 8). The optical flow tracks measure the geodesic cellular surface speed. To measure specifically the component corresponding to mean lateral ruffle travel speed we must project the 3D optical flow velocities along the same plane as the flat cell bottom. The volumetric imaging of cells is acquired on a cover glass tilted at ≈45 0 (Extended Fig. 5c, left). Direct 3D plane-fitting to ( , , ) to determine a precise angle is sensitive to outlier points and shapes deviating from an elongated ellipsoid. Thanks to the ( , ) representation of u-Unwrap3D mapping the curved proximal cell surface and the planar cell bottom to the upper and lower half of the unwrapped 2D image respectively, we could readily annotate in 2D ( , ) the cell bottom and fit a 3D plane to only this surface patch in ( , , ) (Extended Fig. 5c, middle, right). Similarly, we could gate ROI tracks in ( , ) and compute the lateral speed only for those associated the lamellipodia and lamella surface. Doing so we found two populations in the speed histogram (Fig. 5d). Visualizing the speed on ( , , ), the faster population corresponds to higher curvature ruffles with speeds ranging from 2-10 μm/min, and an average speed of 4.2 μm/min (Suppl. Video 8). This is at least two times faster than the slower population surface retrograde actin flow ranging from 0-1 μm/min, which are consistent with the flow speeds we used to measure by 2D fluorescent speckle microscopy in the lamella of epithelial cells 85 . This result suggests ruffles are actively produced and transported across the cell surface. To assess the synchronicity of actin and ruffles we extracted distortion-corrected timeseries of TC and mean curvature, by sampling and averaging the respective values within a spatial window of 23 x 23 centered around the ( , ) coordinates of all ROI tracks on the lamellipodia and lamella surface. Averaging the temporal crosscorrelation curves (mean±95% confidence interval) of individual ROI tracks we find a significant positive instantaneous (lag=0) correlation of 0.2 (Fig. 5e, left). Plotting the instantaneous (lag=0) correlation values of individual ROI tracks as a function of the mean value of the same ROI track and visualizing the instantaneous correlation on ( , , ), we found that ruffles with higher positive surface curvature are more temporally correlated with TC intensity (Fig. 5e, middle, right). Altogether our results show that ruffles are highly dynamic, transient protrusions that actively migrate on the cell surface and driven by locally enriched filamentous actin. Our u-Unwrap3D framework provides now the platform for systematic investigation of the mechanisms that drive and regulate these dynamics.

Discussion
Analyzing the spatiotemporal organization of molecular distributions and signaling activities on cell surfaces in 3D has been limited by the lack of methods to represent, track and process these dynamics. Here we introduce a surface-guided computing framework, referred to as u-Unwrap3D, to bijectively map a genus-X Cartesian 3D surface to equivalent surface and volume representations that are optimally suited for a distinct analytical task. The mappings rely on two critical insights: i) the engineered surface deformation of ( , , ) to generate a genus-0 ref ( , , ) for which a 3D spherical parameterization exists; ii) a novel, efficient algorithm to relax geometric distortion on the 3D sphere in a bijective and tunable manner. Insight i) is fundamental to allowing ref ( , , ) to serve as a representative proxy that captures all the salient surface features of the genus-X ( , , ) surface and insight ii) to preserving this property when ref ( , , ) is unwrapped to the 2D plane. Using diverse cell examples we validated that u-Unwrap3D could be widely applied (>90% of cases) in a manner robust to the input surface mesh quality and that it accurately captures the cell geometry to transfer salient surface features, morphological or molecular, between all representations. We note that this 90% is an underestimation of the applicability. The validation dataset was assembled to be deliberately heterogeneous and was not segmented with the downstream aim of surface mapping. In practice, segmentation algorithms continue to improve and postprocessing techniques can be used to create improved surface meshes. Moreover we can leverage surface meshing algorithms that are more sophisticated than marching cubes, such as dual contouring 86 and shrink-wrapping 87 , to guarantee watertight mesh creation. For timelapse acquisitions, as we showed in Fig. 4,5, the situation is even simpler. Only one timepoint or average surface is needed to generate a single common ̅ ref ( , , ) to unwrap all timepoints. u-Unwrap3D puts in place a generic platform for the spatiotemporal processing of unconstrained cell geometries. u-Unwrap3D is applicable for arbitrary genus-X 3D surfaces, ( , , ) wherever a genus-0 reference surface, Whilst we can increase the range of morphological hole closing in the latter case, we restrict dilation to a maximum 3-5 voxels to minimise smoothing out protrusion features in ref ( , , ). Mesh surgery methods have been developed in computer graphics to make non-watertight 3D meshes watertight 88,89 , and fix imperfections like holes and handles to reduce genus 90,91 to generate higher quality 3D meshes. These have yet to be applied and fully tested for complex cell surfaces. Future development will investigate how to apply such procedures to allow u-Unwrap3D to be applied to input meshes of any quality. In our data we have paid attention during the acquisition to generating sufficient foreground-to-background contrast for reliable surface segmentation, thus minimizing mesh defects.
The concept of geometrical reduction of 3D into 2D geometry through the choice of an optimal coordinate transformation has long existed in mathematics and physics to simplify mathematical manipulation and plotting. For example, parametric coordinates describe the sphere, the cylinder, Mobius strip and helicoids amongst others 92 . In computer graphics this is realised in the common practice of mapping between surfaces through simpler intermediaries: for example, the texture mapping of arbitrary surfaces by optimal surface cutting and mapping of the cuts into individual 2D shapes 57 and ( , ) surface parameterization by cutting and gluing individually mapped 2D planar patches 93 or mapping to canonical shapes such as the triangle 94 , plane 39 or polyhedra 95 to minimize distortion. u-Unwrap3D draws inspiration from this thinking. Through the availability and development of rationalized multiple 3D-to-2D representations, u-Unwrap3D projects analyses that would otherwise require specialised mathematical operations into a sequence of simpler, computationally tractable procedures for 3D mesh processing, image processing and machine learning with specific consideration for single cell biology. Unlike computer graphics benchmarks, surface protrusions are irregular, high curvature and dense. First, we specifically chose representations that map the whole cell surface with well-behaved topologies such as the sphere and plane and designed a relaxation scheme to guarantee interpolation between the minimal conformal and area distortion. This bypasses the numerical instabilities of stitching multiple surface maps and only quasi-conformal mappings for the majority of literature methods. Second, although ref ( , , ) was conceived as a mathematical trick to enable spherical parameterization, because it is derived explicitly from ( , , ) and not as a canonical shape, it serves biologically to decompose ( , , ) as the 'sum' of a smooth cell cortex and surface protrusions. Our results suggest the introduced ref ( , , ) and topo ( , , ) could open up new opportunities to model and quantify the interplay of dynamic membrane morphology and associated signals with more volumetric nuclear/cytoplasmic signalling. Lastly, we designed mappings to underlying representations, including Cartesian 3D, 3D sphere and 2D plane that are standard inputs in computer vision and machine learning.
With u-Unwrap3D standard computer vision, machine learning methods become directly applicable to computing tasks on rugged complex surfaces. More recently, research into combining different geometric representations are state-of-the-art in addressing complex computational problems such as multiple 2D image views to inform 3D mesh reconstruction 96,97 , or 3D mesh vertex coordinates with 2D unwrapped images for feature extraction 98,99 . u-Unwrap3D is fully complementary to these research developments -unifying the different representations into a single surface-guided computing framework for downstream analysis. The u-Unwrap3D framework is made available as a Python library. The resources and validation provided by this work will aid the cell biology community to generate testable hypotheses of the spatiotemporal organization and regulation of subcellular geometry and molecular activity.

Acknowledgments
We thank Hanieh Mazloom Farsibaf and Zach Marin for discussions on mesh processing. We also thank Qiongjing Zou for putting the Python library into a continuous integration coding framework. Funding for this work in the Danuser lab was provided by the grant R35 GM136428 (NIH). AW is a fellow of the Jane Coffin Childs Memorial Fund. GMG is a Damon Runyon Fellow (DRG 2422-21).

Author Contributions
FYZ conceived and developed u-Unwrap3D and conducted the analyses. MD generated the uShape3D protrusion instance segmentations and surfaces used in validation. AW acquired timelapse of the blebbing MV3 cell. GMG, BJC and BC acquired timelapse of the ruffling SU.86.86 cell. FYZ and GD wrote the manuscript with input from all authors. GD provided supervision and obtained funding.

Competing Interests
The authors declare no competing interests.

Data Availability
All data presented herein are available from the corresponding author upon request.
Code Availability u-Unwrap3D is open source, developed as a Python library and is available at: https://github.com/DanuserLab/u-unwrap3D.

u-Unwrap3D framework
Following we describe the algorithms underpinning each step of u-Unwrap3D depicted in Fig.1.
Step 1: ( , , ) to ref ( , , ) Conformalized mean curvature flow (cMCF). cMCF 49 modifies the mean curvature flow (MCF) to avoid the formation of pinches and collapsed vertices that compromise bijectivity and cause early flow termination for watertight meshes with high curvature features. We find cMCF also reduces the size of holes and handles in . The integrals, = ∫ ⋅ and = ∫ t (∇ t , ∇ t ) are called the mass ( ) and Laplacian ( ) matrices at time, .
Substituting this notation and rearranging, the linear algebra MCF equation is is then computed from ( ) by direct matrix inversion. cMCF modifies the MCF equation above by using the Laplacian matrix at time = 0, 0 for all timepoints. The Laplacian matrix is a measure of stiffness between local mesh faces, see active contour cMCF below. Using 0 for all timepoints instead of recomputing implicitly constrains mesh faces to retain the same aspect ratio and this conformalizes the flow.
We use the libigl library 101 with the cotangent Laplacian and barycentric mass matrix as the default implementations of , respectively. We improve the numerics of solving cMCF by normalization of the surface area and recentering of vertex coordinates at the origin after each iteration as recommended in Alec Jacobson's blog post (https://www.alecjacobson.com/weblog/?tag=mean-curvature-flow). u-Unwrap3D implements for optional usage the robust Laplacian of Sharp et al. 102 instead of the cotangent Laplacian, which can improve numerical stability.
Automatic stopping criterion for cMCF. The ideal reference surface ref ( , , ) for topographic representation is the cortical cell shape without protrusions. We find that this corresponds to finding the 'elbow point' in the mean absolute Gaussian curvature (Fig. 1a) and not the convergence limit of cMCF which is the sphere 49 . The difference in the mean absolute Gaussian curvature over vertices, Δ| | ̅̅̅̅̅ = | | ̅̅̅̅̅ − | −1 | ̅̅̅̅̅̅̅̅ between successive iterations − 1 and is used as an automatic stopping criterion for cMCF, stop = max ( min , ). min is a user-specified minimum iteration number and is the first iteration for which Δ| | ̅̅̅̅̅ exceeds a userspecified threshold, Δ| | ̅̅̅̅̅ > Δ ℎ ℎ . We compute the discrete Gaussian curvature, at a vertex , given by the vertex's angular deficit 103 , = 2 − ∑ ∈ ( ) where ( ) are the triangles incident on and is the angle at vertex in triangle .
Mesh voxelization and remeshing. Mesh voxelization converts a surface mesh ( , , ) to a binary volume image where individual voxels are either 1 if they are interior to the surface or 0 if exterior. To do this we create a × × voxel volume image larger than the surface with voxels initialised to 0. We then set the intensities of all voxels indexed by the mesh ( , , ) barycenters to 1, i.e. ( , , ) = 1. To ensure a closed binary volume with all interior voxel intensities = 1, ( , , ) was iteratively subdivided by replacing each triangle face by four new faces formed from adding new vertices at the midpoint of every edge until the mean triangle edge length is < 1 voxel. We use the barycenter coordinates of the final mesh final ( , , ) to set the binary voxel values. In case of small holes in ( , , ) that would prevent a closed binary volume by binary filling only, was first dilated with a ball kernel, then binary infilled, and lastly binary eroded with a ball kernel. The postprocessed was meshed using marching cubes 104 followed by construction of an approximated centroidal voronoi diagram (ACVD) 105 to produce a lower genus remeshing of ( , , ) with approximately equilateral triangle faces. If the input ( , , ) is a smooth shape with only small holes or handles such as the cMCF intermediary ref ( , , ) the proposed voxelization and remesh yields a genus-0 mesh.
Step 2: ref ( , , ) to 2 ( , , ) Quasi-conformal spherical parametrization of genus-0 closed surfaces. For genus-0 closed surfaces the uniformization theorem 52 guarantees the existence of a conformal map onto the unit sphere, 2 . For a closed orientable surface such as ref ( , , ) we can compute the genus, from the Euler characteristic, = 2 − 2 = # − # + # . If = 0 we applied the method of Choi et al. 39,51 which uses the theory of quasiconformal composition to ensure a bijective spherical parametrization with bounded conformal error (i.e. quasi-conformal). In practice, we found conformal errors = 0 (Extended Fig. 3).
Step 3: 2 ( , , ) to Ω 2 ( , , ) Equiareal spherical parameterization by mesh relaxation. We iteratively advect the vertex coordinates of 2 ( , , ), whilst preserving face connectivity and the spherical shape to minimise the per face area distortion factor, . The magnitude and direction to advect each vertex, the vector field ⃗ was found as the solution to the linear heat equation 27  This mesh relaxation bijectively diffuses the area distortion scalar factor on the sphere surface in a stable manner until a triangle face collapses, that is when an interior angle = 0. The extent of area relaxation is determined by the input mesh quality. We find that if ref ( , , ) is regular, with near-equilateral faces and has a minimal number of vertices (> 40k for a pixel resolution of 0.104μm), 2 ( , , ) could be stably relaxed to an equiareal spherical parameterization, Ω 2 ( , , ), without triangle collapse.
Step 4: Ω 2 ( , , ) to ( , ) Automatic determination of unwrapping axis using weighted principal components analysis (PCA). Mapping surface features of interest with minimal geometrical distortion into ( , ) is equivalent to finding an optimal north-south unwrapping axis for Ω 2 ( , , ). This optimization is solved by weighted PCA. Let = ( , , ) denote the coordinate of vertex on Ω 2 ( , , ) and the vertex weight, a score of the importance of mapping this vertex with minimum geometrical distortion. The 3x3 weighted covariance matrix, = ( )( ) over all vertices captures the spread of the weight over the sphere. Eigendecomposition applied to the symmetric matrix finds the principal orthogonal directions of variance given by eigenvalues = [ 1 , 2 , 3 ], 1 ≥ 2 ≥ 3 and eigenvectors = [ , , ]. The eigenvalue captures the concentration of the weight in the direction of the corresponding eigenvector. The optimal north-south unwrapping axis is the smallest eigenvector, . To unwrap with respect to we rotate the vertex coordinates so that is the new z-axis. As the eigenvector matrix is orthonormal and thus a 3D rotation matrix, is the rotation matrix that maps the x-axis, ( UV-mapping the unit sphere. We construct an equidistant UV unwrap of the unit sphere where , the column coordinate equidistantly samples the circumference of the sphere, a total length 2 and , the row coordinate equidistantly samples the arc from north to south pole, a total length . This specifies a × 2 pixel UV image with as a user-defined size. By default = 256 pixels. The UV mapping is constructed by pullback. . Note, setting , , to be the vertex coordinates of any mesh bijective to Ω 2 ( , , ) e.g. ( , , ) in the direct unwrapping case produces similarly the corresponding uv-coordinate mapping for that mesh. The weights , , is also used to map any other vertex associated quantities, i ( Ω 2 ( , , )) such as curvature to i ( ( , )), with vector-valued vertices , , replaced now by scalar values. The construction of the UV map as described above replicates the first and last column of the resulting UV image. For applications such as texture mapping and active contour cMCF (see analysis of ruffles in Datasets section) where the UV image grid needs to be converted to a triangular mesh and the image boundaries 'stitched' together we instead use a × 2 + 1 pixel UV image.
Mesh displacement by active contour cMCF Active contours, or 'snakes' 108 , define the boundary of an image region by minimizing its contour energy, . The contour energy is the sum of an internal, and an external energy, , ( , ) = ( , ) + ( ). The internal energy is set by = ∫ | ′ | 2 + | ′′ | 2 , where the number of ′ denotes the order of the spatial derivative. Here, the first term is the tension and the elasticity of the contour. The second term is the stiffness and the rigidity of the contour. The external energy is set by = −∫ , where is an attractor image for the contour. Minimizing is equivalent to solving the Euler-Lagrange equation, ′′ − ′′′′ = −∇ or in matrix form, + ∇ = 0, where prescribes the constant coefficients for computing the second and fourth order derivatives by finite differences. Given a vertex position ( ), the next position, ( + 1) is computed that minimises the residual error + ∇ using gradient descent and backwards Euler is ( + 1) = ( ) − ( ( + 1) + ∇ ) and the linear system is ( + ) ( + 1) = ( ) + ∇ , where is the identity matrix. If = 0, only comprises the second order coefficients associated with ′′ , and we have ( + −1 ) ( + 1) = ( ) + ∇ or equivalently, ( − ) ( + 1) = ( ( ) + ∇ ) which is identical to the cMCF equation with the mass matrix, , the Laplacian matrix, = and = in response to an external force, ∇ . This general equation can be solved by direct matrix inversion to move surface meshes diffeomorphically. We refer to this as active contour cMCF in this paper. To evolve mesh vertices, ( ) normal to the surface in equal -steps, we set = Φ, the signed distance function, and solve ( − ) ( + 1) = ( ( ) + ∇Φ) iteratively, with ∇Φ evaluated at ( ) for each iteration. A positive moves a cell surface mesh normally outwards from the cell and a negative moves the mesh normally into the cell.
Quantification of geometric deformation errors for meshes Surface mappings do not conserve local geometrical measures like angles, edge lengths and face area. Quantification of the distortion in these measures enables a task-specific optimization of the mapping and correction of statistical measurements made on the mapped surface. There are two primary geometric distortions to quantify; conformal and area distortion error (Extended Fig. 1). An isometric deformation is one with no error; both conformal and area distortion errors are 0.
Conformal error. The conformal or quasi-conformal error, measures the extent the shape of a mesh element , e.g. a triangle face, is stretched. It is 0 if the relative distances between vertices and the angles between edges are preserved after the mapping. We compute of mapping triangle Δ to Δ in 3D by first isometrically projecting all triangles into 2D. Let = ( 1 , 1 , 1 ), = ( 2 , 2 , 2 ), = ( 3 , 3    . We use this metric with = 1, which measures the uniformity of stretch distortion over a surface. This error yields near-equiareal spherical parametrization (Extended Fig. 3b).
Isometric error. We observed tradeoff of conformal error, and log area distortion, log on a similar magnitude scale such that their summation has a unique global minima (Extended Fig. 1e). We thus define an isometric error metric, (1 − ) 1 2 + ( ) log with a constant ∈ [0,1] to weight the relative importance of jointly minimizing conformal and area distortion error. We use = 0.5 in Extended Fig. 3b.

Assessment of geometrical difference between two meshes
Four metrics were used to assess the difference between two meshes 1 and 2 possessing different number of vertices and faces; chamfer distance (CD), Wasserstein-1 distance ( 1 ), the difference in surface area (Δ ) and the difference in volume (Δ ). CD is the mean Euclidean distance between all vertices of 1 when matched to the nearest vertex of 2 and vice versa, Wasserstein-1 distance ( 1 ) or Earth-mover's distance (EMD) is the minimum total area weighted distance of 1-to-1 matching vertices on 1 and 2 . 1 accounts for the area of triangle faces and is minimal if the vertices of 1 is a uniform sampling of 2 or vice versa. Exact computation of 1 is impractical, even for small meshes. We compute 1 using the sliced-Wasserstein, 1 approximation, which uses random spherical projections to sum multiple 1D EMD distances 64 . Specifically we use the ot.sliced.max_sliced_wasserstein_distance function from the Python POT library with 50 projections and average the result from 10 evaluations to report an estimate. The difference in total surface area is Δ = 1 − 2 and is Δ (%) = 1 − 2 2 when given as a percentage. Total surface area was computed as the sum of individual triangle areas. The difference in total volume is Δ = 1 − 2 and is Δ (%) = 1 − 2 2 when given as a percentage. The volume of a mesh was computed as the number of voxels in its binary voxelization computed as described above. We used the minimal possible dilation ball kernel size to ensure a correct volume computation -visual checking of binary voxelization and value at least 3x surface area. We use 1 = topo ( , , ) and 2 = ( , , ) to compute the metrics of Extended Fig. 2,3.

Reference surface inference for measurement of protrusion height
An optimal reference surface for protrusion segmentation must be a ( , ) parameterized surface i.e.  , ) with identical ( , ) coordinate. For both, the mean height, ℎ ̅ is the threshold to give the initial binary segmentation, ( ( , , )) = ℎ ≥ ℎ ̅ . We postprocess by applying graph connected component analysis to remove small segmented regions with surface area < 200 voxels 2 ; diffusing the result using twoclass labelspreading 110 with an affinity matrix, , for 20 iterations, clamping ratio 0.99, and binarizing the label probability with threshold of 0.25 at the start of each iteration. Lastly, any remaining small regions with surface area < 500 voxels 2 was removed. The affinity matrix 69  Topography guided instance segmentation of protrusions Individual protrusions are segmented by identifying high curvature protrusive features and applying connected components analysis. We compute the topographic mean curvature ( ( , , )) = − 1 2 ∇ ⋅̂ with the normal, ̂ given by the unit gradient of the signed distance transform of the binary topographic volume of the cell, ( ( , , )). We compute a binary subvolume restricted to the surface, ( ( , , )), the intersection of the morphological dilation of ( ( , , )) with ball kernel size 2, and the morphological erosion of ( ( , , )) with ball kernel size 2. To identify high curvature surface regions, ℎ ℎ ( , , ) for lamellipodia, we concatenate ( ( , , )) Gaussian smoothed with = 1,3,5 as a 3-dimensional feature for all voxels in ( ( , , )) and apply Gaussian mixture model (GMM) clustering (# classes = 3), keeping the class with the highest mean . To identify ℎ ℎ ( ( , , )) for blebs and filopodia which are circular and smaller, we use ( ( , , )) Gaussian smoothed with = 1 as a 1-dimensional feature for all voxels in ( ( , , )) and apply kmeans clustering (# classes = 3), keeping the class with the highest mean . For efficiency, both GMM and kmeans clusterers are fitted on a random sampling of 10,000 surface voxels. Small regions with < 500 connected voxels are removed. Connected component analysis labels each disconnected region in ℎ ℎ ( ( , , )) as individual protrusions, ( ( , , )). We expand labels by 3 voxels and transfer the segmentation to the surface mesh, ( , , ) by interpolation at the vertex coordinates, ( ( , , )) for further surface-based processing. We first apply the binary protrusion segmentation above, ( ( , , )) to ( ( , , )), taking the intersection and keeping segmentations with size > 100 voxel 2 Cartesian 3D surface area. We diffuse segmentation labels with labelspreading, clamping ratio 0.99 for 10 iterations, with affinity matrix , = 0.9 as above. We do not rebinarize the label probability at the start of each iteration. Finally, we apply ( ( , , )) to the diffused segmentations to get the final instance segmentation labels, ( ( , , )). Equiareal disk parameterization by mesh relaxation. We relax the conformal disk parametrization whilst preserving the boundary topology using the area-preserving flow method 112 . We solve Poisson's equation to compute the smooth vector field for diffusing the area distortion and explicit Euler integration to advect vertex points iteratively with Delaunay triangle flips. The extent of area relaxation achieved is determined by the mesh quality and number of vertices with respect to the extremity of local area distortion. In general, relaxation was less stable compared to our relaxation for spherical surfaces above. For thin and long protrusions, prior downsampling and uniform remeshing of the protrusion submesh was necessary to enable full area distortion relaxation.

Direct 2D unwrapping of protrusion submeshes
To convert a unit disk parameterization to an x pixel image, we 'square' the disk using the elliptical grid mapping formula 113 , multiply the resulting vertex coordinates by /2 and interpolate the coordinates and associated vertex quantities onto a x pixel integer grid. This gives similar results to but is significantly faster than solving the Beltrami equation 94 .

Refinement of undersegmented blebs
Given ( , , ) and the vertex ids corresponding to protrusion , ( ) we first impute any small holes in the segmentation; inner vertices not assigned to protrusion but should be in order to ensure the protrusion submesh, protrusion ( , , ) is a genus-0 open surface. We do this by applying graph connected component analysis on the submesh formed by all vertex ids not part of protrusion , { }\{ } . Any component with number of vertices <10% the total surface area of ( , , ) is assigned to protrusion to form { } impute . The submesh protrusion ( , , ) is formed from { } impute . We downsample protrusion ( , , ) by ¼ the number of vertex points and remesh using ACVD as described above both for computational efficiency and to get a higher quality mesh, protrusion ( , , ) required for computing the intermediate equiareal disk parameterization for a final square parameterization. protrusion ( , , ) is directly unwrapped to a 2D 128 x 128 pixel square image as described above. Positive curvature 'seed' regions are identified by thresholding the mean curvature mapped to 2D, ( protrusion ( , , )) > ℎ ℎ with a global threshold and then applying morphological closing, disk kernel radius 1 pixel. To classify regions as having negative, flat and positive mean curvature, 3-class Otsu thresholding was applied to ( ( , , )) to give two thresholds. All regions with mean curvature greater than the higher threshold ℎ ℎ were positive curvature. Undersegmented blebs correspond to a binary composed of conjoined pseudo-circular regions. We use the gradient watershed 114,115 on the Euclidean distance transform of the high curvature region binary to automatically separate conjoined blebs without seed markers. Mesh matching and interpolation was used to map and segmentation labels between protrusion ( , , ) and protrusion ( , , ) . The refined segmentation were mapped as seed labels from protrusion ( , , ) for every protrusion back to ( , , ). The revised seed labels were then diffused across ( , , ) using the combined geometrical and convexity affinity matrix from above for 10 iterations with = 0.99 . The binary protrusion segmentation from above is applied, and any segmentation with Cartesian 3D surface area < 10 voxels 2 removed to give the final refined protrusion segmentation instances.
Topography guided decomposition of cell surface The instance protrusion segmentation ( ( , , )) above assigns a unique protrusion label ID to each vertex of the surface mesh, ( , , ). We use this surface-based ( ( , , )) as seed labels to partition the total internal cell volume into the volume space unique to each protrusion , ( , , ) and the reference cortical cell volume, ( , , ). This is done in three parts; the construction of ( , , ), volume propagation of ( ( , , )) , and using the previous two parts to volumize individual protrusions to obtain ( , , ) (Extended Fig. 4e-g). Volume propagation of surface-based instance protrusion segmentation. The surface-based protrusion segmentation, ( ( , , )) is converted to voxel-based by setting the value of the voxels corresponding to the integer discretized ( , , ) coordinates to the matching protrusion label ID. We expand the labels by 3 voxels using the Python Scikit-Image skimage.segmentation.expand_labels function and mask with the topographic binary cell volume ( ( , , )) to get the initial topographic volume protrusion segmentation, ( ( , , )) with only the surface of protrusions labelled. We apply marker watershed segmentation slice-by-slice to propagate labels laterally into the protrusion volume within a slice and labels from previous slices, from the top, = + to the bottom, = − of the topographic volume. At a slice = 0 , we use the Euclidean distance transform of ( ( = 0 , , )) for watershed with the seed markers given by the labels of the previous slice, = 0 + 1 combined with the current labels of ( ( = 0 , , )) at slice = 0 . In combining labels, the labels of the previous slice = 0 + 1 takes precedence and overwrites the label of ( ( = 0 , , )). The result, ( ( , , )) assigns a protrusion ID to all voxels in the entire cell volume ( ( , , )) (Extended Fig. 4f). where on the boundary, we use the mass, boundary and Laplacian, 0 boundary matrix defined for a 2D line and mesh , 0 mesh is the mass and Laplacian matrices defined for a 3D triangle mesh.
[ | ] is used to denote the augmented matrix formed by appending the columns of matrix and . We solve for the vertex position at the next timepoint ( + ) as with cMCF by direct matrix inversion.

Surface curvature measurement
The mean curvature, of a 3D surface was measured as the divergence of ̂, the unit surface normal 103 , = − 1 2 ∇ ⋅̂. The surface mesh is voxelized to a binary volume, and ̂ is computed as the gradient of the signed distance transform of with the Euclidean distance metric. computed in this manner as opposed to from the mesh directly using discrete differential geometry 103 or quadric plane fitting 117 which is less affected by the number of mesh vertices or the mesh quality.

Mesh quality measurement
The radius ratio = 2 , defined as twice the ratio between inradius and circumradius was used to measure the face quality for a triangle mesh in Extended. Fig. 2,3. It is a mesh quality measure in the sense that the radius ratio obtains its maximum value of 1 for an equilateral triangle; the shape which jointly maximizes all internal angles and gives the best conditioning number for the mesh Laplacian matrix 118 .

Surface rendering
Triangle meshes were exported from Python using the Python Trimesh library into .obj mesh files and visualized in MeshLab 119 . Volumetric images were rendered in Fiji ImageJ through the volume viewer plugin, and intensities were contrast enhanced for inclusion in the figures using Microsoft PowerPoint. The local surface maximum intensity projection image of Fig. 4c was produced by extending z-axis (depth) rays at every xy pixel, and taking the maximum intensity of voxels within ±9 voxels (±1μm) of the cell surface.

Cell morphology validation dataset
To validate u-Unwrap3D (Fig. 2,3 and Extended Fig. 2-4), we used 66 cell surfaces segmented and surface protrusions classified using u-Shape3D and acquired from high resolution light sheet microscopy 4,120 as previously described 45  Step 5, an upsampling factor of 3 for binary voxelization, = 0.5 voxel steps, a specified = 40 steps and 2D robust smoothing 121 with smoothing factor = 50 for each iteration; for Step 6, binarization of the topographic 3D mapped binary cell segmentation with a threshold of 0.5, then Gaussian smoothing = 1 and initial marching cubes meshing at isovalue 0.5, ACVD remeshing with number of clusters = 50% the number of vertices in the marching cubes mesh.

3D timelapse lightsheet imaging and analysis of blebs
Cell culture and timelapse imaging. All details of the cell line creation, culture and imaging of the MV3 GFPexpressing melanoma cell movie in Fig. 4 were previously published 38 . The movie is a total of 200 frames acquired at a frequency of 1.21 s per frame. Each frame is a 104 x 512 x 512 size 3D volume with a voxel resolution of 0.300 x 0.104 x 0.104 μm.
Cell segmentation and surface meshing. The 200 timepoints were spatiotemporally registered volumetrically to the first timepoint, = 0 as previously described 38 . The cell surface at = 0 was segmented using a multilevel method 122 involving local contrast enhancement, deconvolution and edge enhancement and surface meshed as described above to obtain =0 ( , , ) . Images were deconvolved using the Wiener-Hunt deconvolution approach 123 with our previously published point-spread function 45 .The surface mesh at all subsequent timepoints, ( , , ) were reconstructed using the non-rigid registration deformation field from volumetric registration 38 . The vertex Septin intensity was calculated by extending from the surface a trajectory to an absolute depth of 1 μm along the steepest gradient of the distance transform to the mesh surface, and assigning the 95 th percentile of intensity sampled along that trajectory to the originating vertex to capture the systematically brightest accumulation of Septin signal in the cortical shell. The raw Septin intensity suffers decay from bleaching. We simultaneously normalized and corrected the vertex Septin intensity by computing a normalized Septin intensity as the raw intensity divided by the mean Septin intensity in the whole cell volume at each timepoint.
u-Unwrap3D analysis. We computed a mean surface mesh, ̅ ( , , ) from all ( , , ) as the input surface to u-Unwrap3D. This was done by surface meshing the mean binary volume over all binary voxelizations of individual ( , , ) at an isovalue of 0.5. u-Unwrap3D was applied to ̅ { ( , , ) to create a common static ( , , ) coordinate space that all ( , , ) is mapped to in Step 5 of u-Unwrap3D to generate ( , , ). u-Unwrap3D was run with the same parameters for all steps as for blebs in the validation dataset, except for the following modifications: step 1, the same automatic stopping iteration number but +5 steps, and ACVD with 10% of the marching cubes mesh to get a smoother ref ( , , ); step 4, a smaller 512 x 256 pixel size ( , ) grid and not using the unwrapping axis inferred by curvature-weighted PCA -this axis passed through a large bleb and affected tracking; step 5, = 96 steps -a total of 5μm. Topographic cMCF with the robust mesh Laplacian 102 , mollify factor = 1 × 10 −5 , = 5 × 10 4 was applied to each ( , , ) for 10 iterations to compute the corresponding ( , ).
Bleb segmentation and tracking. Blebs were segmented from ( , , ) at every timepoint using the instance segmentation algorithm with refinement for undersegmented blebs as described above. In computing the binary protrusion segmentation we use a downsampling factor of 4 due to the smaller 512 x 256 pixel ( , ) grid and diffuse the segmentation for 5 iterations as the blebs were smaller than the validation dataset. The segmented 512 x 256 ( , ) bleb images, ( ( , )) were padded 50 pixels on all four sides respecting spherical topology. This is done by periodic padding along the -axis. For the -axis, we pad the top of the image by reflecting the pixels with respect to the first image row (i.e. all pixels in row 2 to row 51) and then flipping in the -axis. Similarly, the bottom is padded by refecting the pixels with respect to the last image row (i.e. all pixels in row 2 to row 511) and then flipping in the -axis. For each unique bleb in every timepoint, we computed the bounding box of the bleb given by top left, ( , ) and bottom right ( , ) coordinates. The bleb bounding boxes were tracked using an optical flow assisted bounding box tracker 77 . Boxes were linked over time into tracks using bipartite matching and the intersection over union (IoU>0.25 for valid match) of bounding boxes as the distance function between pairs. To handle large changes in box size, the matching between the current and next frame was carried out on the predicted bounding box coordinates by local optical flow 124 . Optical flow was computed using the mean curvature, ( ( , )) after rescaling ( ( , )) to be an 8-bit grayscale image using the global minimum and maximum curvature values over time. In case of temporary occlusion or missed segmentation, any non-matched blebs were propagated for up to 5 frames (6s) using the estimated optical flow before track termination. Tracks with > 5 frames (6s) and a mean positive curvature, > 0.1 μm −1 were retained as bleb tracks. The coordinates of retained tracks was corrected to account for the initial padding of 50 pixels. To remove erroneous and duplicated tracks, we uniquely match every segmented bleb in each timepoint to a track by IoU. For each track, we then computed the fraction of its lifetime that could be matched to a bleb and removed all tracks for which this proportion was < 50%. Lastly for each track we checked for sudden changes in the bounding box area, which was indicative of an erroneous bounding box in need of substitution by an inferred corrected bounding box. We applied this procedure to each track in order to construct the timeseries of the bounding box area over the track lifetime and compute a smooth reference timeseries using the central moving average with a window of 3 frames. The bounding box at a timepoint is erroneous if the instantaneous difference between the raw and smooth bounding box area > 500 pixel 2 (the mean ( , ) bleb box area is 361 pixel 2 ). The coordinates of a corrected bounding box is inferred from non-erroneous bounding boxes by interpolation using a linear spline. The tracks that remained fully in-focus over its lifetime were retained for analysis.

Bleb timeseries extraction.
We detected ( , ) pixels on blebs by labelling spatially contiguous areas of positive mean curvature based on 3-class Otsu thresholding defining positive, flat or negative curvature. The largest connected component within a bleb bounding box was defined as on-bleb and the remainder area within the bounding box as offbleb. We extracted distortion-corrected average timeseries of bleb area, mean curvature and septin intensity, that is of a scalar quantity, by observing that the mean of over a Cartesian 3D surface area is equivalent to computing a weighted mean over the equivalent ( , ) area, ∬ ( ( , , )) ∬ = ∬ ( ( , )) ∬ . The weight, is the magnitude of the differential area element, = | × | described above.

Bleb event alignment.
Individual blebbing events were detected within a track by applying peak finding after central moving averaging of bleb area timeseries with a window of 3 timepoints. A peak was defined as having a prominence > 0.5 and separated from a neighboring peak by at least 3 timepoints. Individual bleb event timeseries were constructed and temporally aligned using the detected timepoint of maximal bleb area as timepoint 0 and taking a window of 14 timepoints on either side (a total 29 timepoints, 35 s).

3D timelapse lightsheet imaging and analysis of ruffles
Cell culture and timelapse lightsheet imaging. SU.86.86 cells were purchased from American Type Culture Collection (CRL-1837). The cells were transfected with integrating lentiviral plasmids carrying genes for myristoylated CyOFP1 and Tractin-mEmerald. The cells were cultured in RPMI medium supplemented with 10% fetal bovine serum and 1% anti-anti (Gibco 15240062), at 37°C in a humidified incubator and 5% CO2. SU.86.86 cells were imaged on fibronectin-coated coverslips on a custom axially swept light sheet microscope 125 . The microscope detection system comprises a 25X NA1.1 water immersion objective (Nikon, CFI75 Apo, MRD77220) and a 500mm tube lens. The illumination was done through a 28.6X NA0.66 water immersion objective (Special Optics, 54-10-7). The movie analysed in Fig.5 is a total of 30 frames acquired at a frequency of 2.27 s per frame. Each frame is a two-channel 151 x 1024 x 1024 size 3D volume with a voxel resolution of 0.300 x 0.104 x 0.104 μm.
Cell segmentation and surface meshing. All timepoints were rigid registered volumetrically to the first timepoint, = 0 to compensate for drift. The CyOFP1 image was also rigid registered to the Tractin-mEmerald in each timepoint. For every timepoint, the volumetric image was segmented using a multi-level method 122 involving local contrast enhancement, deconvolution and edge enhancement and surface meshed as described above to obtain the surface mesh, ( , , ) . The vertex Tractin-mEmerald and CyOFP1 intensities were calculated by extending a trajectory to an absolute depth of 1 μm along the steepest gradient of the distance transform to the mesh surface, and taking the mean intensity along the trajectory.
u-Unwrap3D analysis. The first timepoint surface mesh was used as the input ( , , ) to u-Unwrap3D to create a common static ( , , ) coordinate space that the surface meshes from all timepoints is mapped to in Step 5 of u-Unwrap3D to generate ( , , ). We use Unwrap-3D with the following parameters for each step: Step 1, cMCF with maximum iterations = 50, = 1 × 10 −5 , automatic stopping threshold, Δ ℎ ℎ = 5 × 10 −5 , ref ( , , ) mesh voxelization with morphological dilation and erosion with ball kernel radius 5 voxels, Gaussian smoothing = 1 of the binary volume and initial marching cubes meshing at isovalue 0.5, ACVD remeshing with number of clusters = 10% the number of vertices in the marchin cubes mesh, and further volume constrained Laplacian mesh smoothing 126 with implicit time integration, time step 0.5 for 15 iterations; for Step 3 area distortion relaxation with maximum iterations = 100, = 0.1, stepsize = 1; for Step 4 we use the binary positive curvature region of ref ( , , ) given by 3-class Otsu thresholding as the weight for determining the unwrapping axis, and use a 1025 x 512 pixel ( , ) grid; for Step 5, for outwards propagation, an upsampling factor of 3 for binary voxelization, = minimum of 0.5 and 1 2 (Δ + Δ ) voxel -steps where Δ , Δ is the mean Cartesian 3D distance of traversing one pixel in , directions and a separable 1D uniform box filter smoother with a window 5 pixels; for inwards propagation, we use active contour cMCF with = 5 × 10 −4 , and robust mesh Laplacian 102 , mollify factor 1 × 10 −5 for better numerical stability with the ( , ) parameterized ref ( , , ) converted into a triangle mesh by triangulating the quadrilateral pixel connectivity and inserting additional triangles to 'stitch' the image boundaries into a spherical topology, (note the latter stitching requires an even number of columns after discounting that the last column is the same as the first column, hence a 1025 x 512 grid) and = (maximum internal distance transform value ) / steps; for Step 6, marching cubes meshing of the topographic 3D mapped binary cell segmentation at isovalue 0.5 following Gaussian smoothing = 1, ACVD remesh with number of clusters = 50% the number of vertices in the marching cubes mesh and retaining the largest connected component mesh. Topographic cMCF with robust mesh Laplacian 102 , mollify factor = 1 × 10 −5 , = 5 × 10 4 was applied to each ( , , ) for 10 iterations to compute the corresponding ( , ).
Optical flow ROI tracking. We used motion sensing superpixels (MOSES) 83,84 in dense tracking mode which automatically monitors the spatial coverage of ROIs and introduces new ROIs dynamically to ensure uniform spatial tracking at every timepoint. We partitioned the image with an initial user-specified 1000 nonoverlapping rectangular regions-of-interest (ROI). Each ROI was tracked over time by subsequently updating its centroid by the median optical flow 124 within the ROI. Optical flow was computed from the Traction/CyOFP1 (TC) signal after rescaling TC( ( , )) to be an 8-bit grayscale image using the video minimum and maximum TC values.
ROI timeseries extraction and cross-correlation. Distortion-corrected average timeseries of mean Tractin/CyOFP1 (TC) and were computed for a track using the same weighted mean as for blebs. A square bounding box of the mean MOSES ROI width centered at the track ( , ) coordinate was used to sample the scalar values at each timepoint. The distortion-corrected timeseries can be treated as standard 1D timeseries. The 1D normalized cross-correlation was thus computed between the distortion-corrected TC and timeseries for individual tracks without modification. ROI cross-correlation curves were averaged at all timelags to derive the mean and 95% confidence interval ROI cross-correlation curve. A deviation of the curve greater than the 95% confidence interval at a time lag of 0 indicated significant instantaneous correlation.
Retrograde actin flow and mean ruffle travel speed. Computing the speed histogram with 25 bins and speed range 0-10 μm/min showed a slow and fast population (Fig. 5d). We inferred the mean speed of the two populations as the two thresholds generated by 3-class Otsu thresholding. The lower and faster of the thresholds are the mean speed of retrograde actin flow and ruffles respectively.
Cross-correlation and curvature relationship. We computed the continuous relationship of mean curvature, and the lag 0 cross-correlation of TC and (Fig. 5e) over ROI tracks using kernel density. Gaussian kernel density with a bandwidth set by Scott's rule was used to derive the joint density distribution of and cross-  Overview of the surface-guided computing framework u-Unwrap3D. a) Overview of the 6 key steps to map an input genus-X Cartesian 3D surface, ( , , ) such as that obtained from surface meshing the binary segmentation of an input 3D volume image, , and associated scalar measurements, ( ( , , )), via a smooth genus-0 reference surface, ref ( , , ), into any of three additional representations; topographic 3D surface, ( , , ), 3D sphere, 2 ( , , ) and 2D plane, ( , ). denotes mesh vertex coordinates, is the mesh mass matrix, is the mesh Laplacian matrix, the step size of area-distortion relaxation, the step size (in pixels) of the propagation distance, the topographic depth (in pixels) and the iteration number. c) Steps to directly remap input genus-0 surfaces without need for a reference surface. In the figure, (⋅), (⋅) denote surface and volume geometries, respectively, in either Cartesian 3D, topographic 3D, or radiusstandardized 3D spherical coordinates; ( (⋅)) and ( (⋅)) denote surface or volumetric signals as a function of a particular surface or volume geometry. and denote mean and Gaussian curvatures respectively. Figure 2. u-Unwrap3D generates a spectrum of equivalent data representations for surface-guided computing. a) Summary of the bijective mappings between the 5 equivalent surface mesh representations generated by u-Unwrap3D. Black bidirectional arrows indicate the mapping algorithms between representations discussed in the text. Grey arrow indicates the direct spherical mapping applicable when the input mesh is genus-0. b) u-Unwrap3D also enables bidirectional mapping of volumetric information between a Cartesian and topographic space relative to a genus-0 reference surface. c) Gallery of equivalent surface representations generated on examples of cell surfaces with blebs, lamellipodia and filopodia. For visualization of the mappings, individual instances of morphological motifs detected by the software uShape3D are color-coded on surface representations. d) Gallery of equivalent volume representations for the same cells shown in c). Volume image intensities were visualized using ImageJ volume viewer and contrast-enhanced to better visualize fine protrusions (see Methods).